Individual and Collective Behavior of Small Vibrating Motors 
Interacting Through a Resonant Plate 
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Abstract 



We report on experiments of many small motors - cell 
phone vibrators - glued to and interacting through a 
resonant plate. We find that individual motors inter- 
acting with the plate demonstrate hysteresis in their 
steady-state frequency due to interactions with plate 
resonances. For multiple motors running simultane- 
ously, the degree of synchronization between motors 
increases when the motors' frequencies are near a res- 
onance of the plate, and the frequency at which the 
motors synchronize shows a history dependence. 

1 Introduction 

Ensembles of oscillators that spontaneously syn- 
chronize have been studied for decades. Biolog- 
ical examples aboundpQ [3J |3], but synchroniza- 
tion occurs in many other systems, including cou- 
pled metronomes [H [5], laser arrays [6|, chemical 
oscillators [7j, arrays of convective cells [8], Josephson 
junctions arrays [5], transport networks|10|. and per- 
haps most notoriously, pedestrians crossing the Mil- 
lennium Bridge in London when it first opened^llj. 
These systems are all examples of populations of sim- 
ilar but not identical oscillators that exhibit the same 
basic patterns of behavior, that (1) they synchro- 
nize spontaneously, without the need for any external 
driving, and (2) as the oscillators' coupling increases, 
their synchronization strengthens. For an overview of 
the topic, see the review by Acebron et al. [12] and 
the popular book Sync, by Strogatz|13). 

The topic of synchronization is much broader 
than the study of many coupled oscillators. In 
an effort to better understand radio tuning, Adler 
studied the synchronization of locking circuits, in 
which a phase-oscillator synchronizes to a periodic 
forcing [14]. Burykin and Buchman discussed the pos- 
sibly lethal outcome of the lack of synchronization 
among organ systems when taking a patient off of 
a mechanical respirator[15i. Gintautus et al. found 
synchronization in so-called mixed-reality states, in 



which virtual and real systems are coupled and inter- 
act in real time|16|. All of these systems exhibit syn- 
chronization in some sense. Although we find these 
systems to be interesting, the work presented here is 
motivated by the many examples listed in the first 
paragraph: spontaneous collective behavior of many 
coupled oscillators in the absence of external forcing. 
In this paper we present yet another system that 
exhibits synchronization: small mechanical vibrators 
coupled through a resonant plate. In addition to 
being inexpensive and easy to study, this system 
provides a unique twist to the standard coupled- 
oscillator problem in that the coupling between the 
oscillators depends on frequency and exhibits a sim- 
ple resonance structure. How does frequency depen- 
dent coupling effect the dynamics of coupled oscilla- 
tors? Unlike most other globally coupled oscillator 
systems, we find history-dependent behavior: an en- 
semble of oscillators shows hysteresis in the frequency 
and degree of synchronization, and an individual os- 
cillator shows hysteresis in its steady-state frequency. 



2 Experimental Setup 

In this work, we study 16 small motors with eccen- 
trically massed rotors. The motors (All Electronics 
Corporation, catalog number DCM-204j) are small 
DC motors, the sort used as vibrators in mobile 
phones. Each motor has a mass oi 3g and is 2 cm 
long. Each motor's rotor has a center of mass that is 
offset from the axis of rotation, with a first moment 
of 0.74 g-mm. Vibrations arise from the rotation of 
this off-center mass. 

To cause the motors to interact, we attach them 
to a mechanically compliant and resonant aluminum 
plate held by clamps as shown in figure [TJ The plate 
is L = 115 cm long, b = 15 cm wide, and 5 mm thick. 
We adjust the linear response of the system by ad- 



^This item is no longer available in the catalog, but simi- 
lar motors can be found in their catalog searching for "motor 
vibrator." 
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Figure 1: A photo and diagram of the experiniental 
setup, as described in the second section. 



Figure 2: Typical time series of a single motor on 
the plate for different voltages. From top to bottom, 
the data correspond to driving voltages are 0.65 T^, 
0.84 y, and 1.05 F. 



justing the location where the clamps hold the plate, 
parametrized by the length a. Although we consid- 
ered various clamp positions, all of the results re- 
ported here are based on a length of a = 12.5 cto, 
for which the system has resonances at 68 Hz and 
100 Hz. 

We measure the plate's vertical acceleration a (t) 
using an accelerometer attached to the plate, a PCB 
353B33. In the diagram shown in figure [U the ac- 
celerometer is depicted by the canister underneath 
the motors. A few typical time series of acceleration 
data due to a single motor are shown in figure [5] The 
sampling rate for these and all other data we dis- 
cus is r = 1000 iJz. The plate is a linear medium, 
so we attribute any observed vibrations either to the 
motors or to background sources, such as building vi- 
brations. We took measurements with a stroboscope 
to confirm the primary frequencies obtained from the 
accelerometer data. In order to reduce spurious fre- 
quencies from the environment, we place the entire 
setup on a foam pad. Although some background 
noise still perturbs the system, these vibrations do 
not dominate the signal reported by the accelerome- 
ter and have frequencies much lower than the motors' 
primary frequencies. 

When the motors operate near a resonance, the 
stroboscope allows us to observe the plate's mode 
shape. We find that both resonances, near / = 68 Hz 
and 100 Hz, do not have any nodes along the array of 



motors and that the displacements of all the motors 
has about equal magnitude. As such, the coupling 
between the motors has no appreciable position de- 
pendence. 

All of the motors operate from a common power 
supply but small variations in each of the motors lead 
to a distribution of motor frequencies for a given volt- 
age. We do not attempt to characterize the distribu- 
tion of motor speeds in any rigorous way since we 
are only working with 16 motors. Stroboscopic ob- 
servations indicate that the frequency distribution is 
approximately unimodal. 



3 Behavior of a Single Motor 

In order to discuss how multiple motors interact we 
must first understand how a single motor behaves and 
the sort of response it produces. 

Figures [5] and [3] demonstrate typicao single mo- 
tor data at vohages V = 0.65 V", 0.84 y, and 1.05 V^. 
Both figures show stable periodic behavior. The spec- 
tra a (/) in figure [3] were computed from two-second 
long data sets a {t) by short-time Fourier Transform: 



a(/)= / a{t)e 



„' 2 TT / 



dt. 



(1) 



^Although the response as a function of primary frequency 
is typical, the primary frequency as a function of vohage is 
not. This is one of our faster motors. 
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Figure 3: Typical Fourier transforms of a single mo- 
tor on the plate for different voltages. The driving 
voltages are 0.65 y (-), 0.84 y (•••), and 1.05 F (- -). 



as implemented with an FFT. Although the motor's 
velocity can drift under special circumstances, the 
narrow widths of the peaks in figure [3] demonstrate 
that a motor's velocity is relatively stable. We iden- 
tified the primary frequency / of the motor in these 
and many other similar experiments by fitting the 
peak to a Lorentzian. When we plot the primary fre- 
quencies verses the driving voltage we obtain the plot 
shown in figure |4lja), which we will discuss in more 
detail below. We also see that the magnitude of the 
response varies substantially for different primary fre- 
quencies. Although the magnitude of the signal can 
be measured by its maximum value, a more robust 
measure is the RMS of the Fourier transform in the 
vicinity of the peak, Mj^ms- The RMS magnitude 
can be plotted against the voltage V , but it is better 
understood as a function of the primary frequency /, 
as shown in figure [5] We will return to the magni- 
tude measurements shortly, but first we will examine 
how the primary frequency depends on the applied 
voltage. 

Although the motor's primary frequency is rela- 
tively stable when the voltage is fixed, figure ID^a) 
shows that the motor's frequency is not a simple func- 
tion of voltage. We can compute the primary fre- 
quency / and the width of the peak 5f by fitting the 
magnitude of the spectrum a to a Lorentzian, 
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Figure 4: Frequency response of a single motor verses 
voltage, both (a) on a resonant plate, and (b) for 
comparison, on a rigid support. The motor studied 
in figure (a) was different from the motor studied in 
figure (b). 



within the vicinity of of the peak. Shown in fig- 
ure IH^a) are the primary frequencies for two different 
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Figure 5: Accelerometer amplitude as a function of 
motor frequency when being driven by a single motor. 
The peaks near 68 Hz and 97 Hz correspond with 
peaks in the support's Green function at the same 
frequencies. 

sets of consecutive measurements, one in which we 
started at y = 2AV and slowly decreased the volt- 
age to 0.6 V (indicated by triangles pointing down- 
ward), and another in which we started the motors 
at V = 0.6 V and slowly increased the voltage to 
2AV (indicated by the triangles pointing upward). 
Although the two measurements demonstrate rela- 
tively good agreement below V — IV and above 
T^ = 2 y, we see a hysteresis between V = 1 and 
2 V. The upward data gets stuck near a resonance of 
the plate. In contrast, similar data taken from a sep- 
arate motor on a rigid support is shown in fLgure|4l(b) , 
and we see that in the absence of resonances, a mo- 
tor's frequency is nearly linear in the applied voltage. 
The marked difference indicates that the motor in- 
teracts strongly with the resonances of the plate, and 
these interactions lead to the hysteresis observed in 
figure mja) . 

The primary frequency's magnitude M^ms shows 
a strong dependence on the primary frequency /, as 
shown in figure O We compute the RMS magnitude 

by 



M, 



RMS 



f+^f 



\h- 



-A/ 



l«(/)l d/, 



(3) 



where the range A/ is a multiple of the width 5f de- 
termined by the Lorentzian fit, equation [51 The val- 
ues obtained for Mrms are largely independent of the 
choice of A/ so log as A/ > 5f. As shown in figure[5l 
the magnitude of the plate's response to a single mo- 
tor is not monotonic in frequency. We can understand 
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Figure 6: Normalized plot of data shown in figure [3] 
The driving voltages are 0.65 V (-), 0.84 V (•••), and 
1.05 F (--). 



this behavior by noting that the plate has resonances 
near 10 Hz and 100 i?z, so the plate will have larger 
accelerations when driven by a motor near these fre- 
quencies than when the motor's frequency is far from 
the resonances. This data was obtained by powering 
different motors - one at a time - at various voltages 
and taking two-second data sets for each voltage. Al- 
though we could seek a relation between the RMS 
magnitude and the applied voltage, figure [S] indicates 
that the RMS magnitude is a function of primary fre- 
quency. Despite overlaying data from motors at var- 
ious different locations on the plate, the magnitude 
as a function of primary frequency is remarkably con- 
sistent, showing that the coupling between the plate 
and the motors for our geometry does not depend 
substantially on the motor's position. Apart from 
the gaps in the data for frequencies just above the 
two peaks, the magnitudes in this plot are equivalent 
to f^ G, where G represents the passive frequency- 
dependent Green function of the system. Calcula- 
tions not shown here indicate that the gaps are due 
to motor-plate interactions that make those frequen- 
cies unstable. 

All of the discussion of data presented so far has 
focused on single motors. Since we use a single ac- 
celerometer to measure the behavior of multiple mo- 
tors acting simultaneously, and since we wish to know 
when two motors synchronize, we must obtain a rea- 
sonable estimate for the number of motors at a given 
frequency. Such an estimate is not trivial: the res- 
onant response of the plate means that one motor 
turning at 95 Hz will produce a much stronger signal 
than many synchronized motors with a primary fre- 



quency of 78 Hz. Our solution to this problem is to 
use figure [5] as a normalization curve. We sample the 
RMS magnitude uniformly - interpolating where nec- 
essary - to obtain normalization amplitudes M(/). 
We then normalize a raw spectrum such as figure [3] 
by dividing the amplitudes of the original spectrum 
by the normalization amplitudes: 



N{f) 



l^(/)l 
MifY 



(4) 



The result of such a normalization scheme is shown in 
figure [6] for the data presented in figure [S] Except for 
the artifacts at / = 50 Hz and 75 Hz associated with 
the signal at F = 1.05 V, the scheme appears to work 
quite well. Even with the artifacts, single motors 
can be easily distinguished and counted, providing 
us with a decent measure of the number of motors in 
the vicinity of a given frequency. 

4 Many Motors on a Resonant 
Plate 

The essential behavior of multiple motors interacting 
on the plate is given in figure [T] These plots are con- 
secutive minute-long measurements that have been 
Fourier transformed and normalized as discussed in 
the previous section. Instead of plotting individual 
spectra, like those in figure [SI we plot consecutive 
spectra by creating gray-scale columns and laying 
them out sequentially in order of applied voltage. 
The difference between figures [71[a) and[71[b) is that 
in the former we started the system at high voltage 
and stepped the voltage down each consecutive mea- 
surement, whereas in the latter we started the sys- 
tem at low voltage and stepped the voltage up each 
consecutive measurement, in a manner similar to the 
two data sets shown in figure 31 Our motors have 
multiple stable states for a given voltage, but gener- 
ally the motors show partial or full synchronization 
near f = 65Hz (1.21/ < V < lAV), partial syn- 
chronization below f = 85Hz (1.2 F < T^ < 1.7 V"), 
and nearly full synchronization above / = 85 Hz 
{V > 1.7 F). In light of the multiple stable states of 
a single motor near a resonance, the multiple stable 
states of the many- motor system between 1/ = 1.2 1/ 
and 1.71/ is not surprising. 

However, the frequency of the resonance that 
causes the hysteresis is surprising. Although the mag- 
nitude measurements in figure [S] clearly show the res- 
onance near 68 Hz, the individual motor's behavior 
shown in figure |4] indicates that the resonance has no 
noticeable effect on the motor's frequency. The pro- 
nounced effect of the resonance in figures [7{^a) and (b) 
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Figure 7: Behavior of many motors on a plate as 
a function of voltage, (a) Behavior as we decrease 
the voltage starting from an initially high value, (b) 
Behavior as we increase the voltage from an initially 
low value. 
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Figure 8: Normalized spectrograms of typical and 
atypical dynamics of multiple motors on a resonant 
plate, (a) Behavior at 1.49 V. (b) Behavior at 1.06 V. 



suggest that a resonance's effect on a motor's steady- 
state frequency depends on the number of motors 
near the resonance. We assert that when we have 
more motors near a resonance, the resonance's pull is 
stronger. To confirm this assertion, note that in fig- 
ure[71^b), between V = IV and 1.1 y, there is a motor 
turning with increasing frequencies, from / = 65 Hz 
to 70 Hz. At 1.15 V" its frequency is 74 Hz. Before 
proceeding to 1.2 V", we forced the motor back down 
to the ensemble near 62 Hz. Had we continued the 
measurements with that motor left unchecked, as we 
did in other measurements, the synchronization at 
62 Hz would have dispersed at V = 1.45 F or 1.5V. 
This difference in behavior confirms that the stabil- 
ity of a resonance depends on the number of motors 
within a vicinity of the resonant frequency. 

Figure [5] shows the dynamics of the motors at a 
fixed voltage over long periods. The stability of the 14 
motors synchronized at / = 80 Hz is typical for syn- 
chronized motors, while the behavior of the barely- 
visible motor wandering near 65 Hz is typical of un- 
synchronized motors. The behavior captured in fig- 
ure [Sljb) is unusually dynamic for steady-state behav- 
ior: it demonstrates nearly all of the typical transient 
behavior that the motors exhibit at higher voltages. 
The most interesting behavior is the synchronization 
and desynchronization of motors spinning faster than 
/ = 55 Hz. In particular, two groups of motors grad- 
ually merge between f = 50s and 70s, and then 
merge with yet another group to form a very large 
group at 130 s. This group is not stable and some of 
the motors split off just before 200 s before re-merging 
about a minute later (not shown). 



5 Conclusion and Future Work 

We set out to answer this question: How does 
frequency-dependent coupling effect the synchroniza- 
tion dynamics of many coupled oscillators? We found 
that the behavior of individual motors and ensem- 
bles of motors interacting with a resonant plate is 
reproducible but history dependent. However, much 
remains to be known. 

In the work presented here, we studied the behavior 
of one motor on a rigid support and others on a res- 
onant plate. In future experimental work, one could 
measure the same motor's angular velocity when it 
is on a rigid support and on a plate to get a direct 
comparison of frequency as a function of voltage. We 
also considered a specific geometry for the plate and 
motors so that all of the motors interacted with the 
plate in nearly the same way. How would the motors 
behave differently if some of them were placed on the 



nodes for a given resonance? 

We have not presented a model for our system, 
but many of the behaviors exhibited in this paper 
seem theoretically tractable. These include the sta- 
bility criteria for a single motor as it approaches a 
resonance from above or settles on a resonance from 
below, the stability criteria for synchronization near 
a resonance, and the steady-state dynamics of many 
motors at weak coupling. Is the behavior at weak 
coupling chaotic? Is it stochastic? What governs the 
time scales of merging and collapsing groups? 

We look forward to studying all of these fascinating 
behaviors. 
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